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For the models of inflation driven by the potential energy of an inflaton field cj>, the covariant 
Galileon Lagrangian (d(f>) 2 04> generally works to slow down the evolution of the field. On the other 
hand, if the Galileon self-interaction is dominant relative to the standard kinetic term after the end 
of inflation, we show that the oscillation of inflaton tends to be violated during reheating. This is 
typically accompanied by the appearance of the negative propagation speed squared c 2 s of a scalar 
mode, which leads to the instability of small-scale perturbations. For chaotic inflation and natural 
inflation we clarify the parameter space in which the violation of inflaton oscillations does not occur. 
Using the WMAP constraints of the scalar spectral index and the tensor-to-scalar ratio as well, we 
find that the self coupling A of the potential V(<f>) = A^> 4 /4 is constrained to be very much smaller 
than 1 and that the symmetry breaking scale / of natural inflation cannot be less than the reduced 
Planck mass M p \. We also show that, in the presence of other covariant Galileon Lagrangians, there 
are some cases in which the inflaton oscillations are not violated even for the self coupling A of the 
order of 0.1, but still the instability associated with negative c 2 is generally present. 
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I. INTRODUCTION 



The idea of inflation was originally proposed to address a number of cosmological problems plagued in standard 
Big Bang cosmology Moreover inflation provides a causal mechanism for the generation of large-scale density 
perturbations from the quantum fluctuation of a scalar field ( "inflaton" ) Q . The resulting power spectra of scalar and 
tensor perturbations are nearly scale-invariant, whose prediction is consistent with the Cosmic Microwave Background 
(CMB) temperature anisotropies observed by COBE [3| and WMAP 

Most models of inflation are based on a canonical scalar field <f> with a slowly varying potential V(<p) (see [|| for 
reviews). For example, the simple power-law potential V((f>) = X(f> n /n (n and A are positive constants) leads to chaotic 
^ ■ inflation for the field value larger than the reduced Planck mass M p i = 2.435 x 10 18 GeV @. The quartic potential 
V(4>) = \(j} 4 /4: is in tension with the WMAP constraints of the scalar spectral index n s and the tensor-to-scalar ratio r 
Q. Moreover the self coupling is constrained to be A « 10~ 13 from the WMAP normalization, which is much smaller 
than the typical coupling scale appearing in particle physics (e.g., A w 0.1 for the Higgs boson Q). 

There are several different ways to reconcile the quartic potential V(cb) = A</> 4 /4 with observations. One of them is 
to introduce a non- minimal field coupling £Rcf) 2 /2 to the Ricci scalar R 0, [HJ. In the limit £ 3> 1 the tensor-to-scalar 
| ratio can be as small as r ~ 10 -3 with n s m 0.96 [ll|, which is well inside the la observational contour [l2[. Moreover 
the self coupling is of the order of A w 10~ 10 £ 2 for £ ^ 1 from the WMAP normalization. If the field 4> is a Higgs 
boson, however, this model is plagued by the problem of unitary violation around the energy scale of inflation [13j j . 
Moreover the non-minimal coupling £i?</> 2 /2 does not necessarily help other inflaton potentials to be compatible with 
^ ; observations [H [U] . 

The second way is to use a non-minimal field derivative coupling to gravity in the form G^d IJl cj)d l/ (j)/(2M 2 ) [Pol ], 
where G^ v is the Einstein tensor and M is a mass scale, (see also Ref. [l|| for the original work). In the regime 
where the Hubble parameter H is larger than M, the evolution of the field slows down due to a gravitationally 
enhanced friction. In this case the potential V(4>) = A^> 4 /4 is compatible with the WMAP constraints of n s and r 
with A - 5.9 x 10~ 32 (M pl /M) 4 [13 • Moreover the mechanism of slowing down the field ( "slotheon" [ijj ) works for 
general steep potentials. For example this mechanism was applied to the potential VUjA = A 4 [1 + cos(</>/ /)] of natural 
inflation, where A and / are mass parameters (l9j . In conventional natural inflation )2Ci the symmetry breaking scale 
/ needs to be larger than 3.5M p i for the consistency with the WMAP constraints (21|, but in this regime standard 
quantum field theory is unlikely to be trustable [22j . In the presence of the field derivative coupling, natural inflation 
can be compatible with the WMAP bounds even for / smaller than M v \ [13, HH . 

For the potential-driven inflation the field self-interaction of the form (d4>) 2 \34> [24T - [26| also leads to the slow 
evolution of inflaton along the potential [27j (see Refs. [28[ for the kinetically driven case). The field equations of 
motion following from the Lagrangian (9</>) 2 D0 respects the the Galilean symmetry 9 M </> — > d^cj) + & M in the limit of 
Minkowski space-time [24[. In a manifold having integrable (covariantly constant) Killing vectors £ a , such a "Galileon" 
Lagrangian is invariant under the curved-space Galilean transformation <p(x) — > <fi(x) + c + c a J XQ £ a , where c, c a , xq are 
constants and a; is a space-time coordinate [l~8| (whose property also holds for the derivative coupling G^ lh ' dgcbdycjA. 
The presence of such a symmetry has an advantage that the theory can be quantum mechanically under control [29J . 
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For the potential V(4>) = A</> 4 /4 the Galileon term (d<fi) 2 E\(j) not only leads to the suppression of the tensor-to-scalar 
ratio compatible with recent observations, but also it gives rise to the coupling A of the order of 0.1 consistent with 
the WMAP normalization [3, H3, 0I3|. Meanwhile it is not clear whether the presence of such a non- linear field 
self-interaction does not disturb the oscillation of inflaton during reheating. Moreover we need to check whether the 
conditions for the avoidance of gho sts and Laplacian instabilities can be avoided after inflation. Since such conditions 
were recently derived in Refs. |3ll-l34l] for the most general scalar-tensor theories having second-order equations of 
motion [3514371 ] , those results can be applied to potential-driven inflation with the Galileon Lagrangian. 

In this paper we study the dynamics of inflation and the subsequent reheating for the potentials V(<fi) = \(j) n /n 
and V(4>) = A 4 [l + cos(</>//)] in the presence of the Galileon Lagrangian (d(f)) 2 U\(f). If the Galileon self-interaction 
dominates over the standard kinetic term after inflation, the oscillation of inflaton tends to be violated for both 
potentials. This is usually accompanied by a negative propagation speed squared c 2 of the scalar mode, which leads 
to the instability of scalar perturbations on smaller scales. The model parameters of the potentials can be constrained 
to avoid the violation of the inflaton oscillation as well as to match with the observational data. For the quartic 
potential V(4>) = A^ 4 /4, for example, the self coupling A is bounded to be very much smaller than 1. In natural 
inflation we show that it is difficult to realize the regime where the symmetry breaking scale / is smaller than M p \. 
We also study the effect of other covariant Galileon terms [25[ on the dynamics of inflation and reheating for the 
potentials V((f>) = X(f> n /n. It is possible to find some cases in which the self coupling of the potential V(4>) = A(/> 4 /4 
is of the order of 0.1, but the violent instability associated with negative c 2 is usually unavoidable. 

This paper is organized as follows. In Sec. |IT]we present the background and perturbation equations for potential- 
driven inflation in the presence of (generalized) Galileon Lagrangians. The spectra of scalar and tensor perturbations 
are given by using slow-roll parameters. In Sec. lIIII we study the models of chaotic inflation as well as natural inflation 
in the presence of the term (d(j)) 2 n\(j) alone. We clarify the viable parameter space in which the oscillation of inflaton is 
not violated during reheating. We also place observational constraints on the inflaton potentials from the information 
of the scalar spectral index n s and the tensor-to-scalar ratio r. In Sees. IIVI and [V] we provide similar constraints 
on the parameter space of chaotic inflation in the presence of other covariant Galileon terms. Sec. I VII is devoted to 
conclusions. 



II. GENERAL FIELD EQUATIONS FOR THE BACKGROUND AND PERTURBATIONS 

We start with the following action 



S = Id xy/—g 



Ml 



-^R + P{<t>, X) - G 3 (0, X)U<j> + U + d 



(1) 



where g is a determinant of the metric g^ , R is a scalar curvature, and 

U = G 4 (<^, X)R + G 4 , x [(D0) 2 - (V^V^XVW)] , (2) 
C 5 = G 5 (0,X)G^(VV» - ^G 5 . x [(D0) 3 - 3(D0)(V AI V^)(V^V» + 2(VV Q 0)(V Q V /3 0)(V< 9 V^)] . (3) 

Here P and Gi (i = 3,4,5) are functions in terms of </> and X = —g fJ ' ,/ d ll (f>d^cf>/2 with the partial derivatives G^x = 
dGi/dX, and G M „ = R^ — <? M „i?/2 is the Einstein tensor (R^ is the Ricci tensor ). T he action (fl]) corresponds to 
the most general scalar-tensor theories having second-order equations of motion [HI HH, H3] ■ This was first discovered 
by Horndcski in a different form [35| . 

We focus on the models in which inflation is mainly driven by a field potential V(<f>), i.e., 



P(<f> 1 X) = X-V( ( f>). 

For the functions Gi (i = 3, 4, 5) we take 

G a (4>,X) = M<j>)X, G 4 (<j ) ,X) = f 4 (q>)X 2 , Gb(4>,X) = h{4>)x 2 
where fi(4>) depend on <f) alone. The covariant Galileon J2f| corresponds to the choice [38[ 



h 



_£3_ 

M 3 



C4 

M 6 



3c 5 

M9 ' 



(4) 
(5) 

(6) 



where c 3 , C4, C5 are dimensionless constants, and M is a constant having a dimension of mass. We derive the 
background and perturbation equations for the general functions © in order to cover the cases discussed in Refs. (27l . 
|39| (such as fd,{4>) <x 4>). After Sec. Mil we mainly focus on the covariant Galileon. 
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A. Background equations 



On the flat Friedmann-Lemaitre- Robertson- Walker (FLRW) space-time with the scale factor a(t) (where t is cosmic 
time) the background equations for the theories described by the functions Q and (JSJ) are given by jU HH 

E x = 3M^H 2 -X-V- 6Hf 3 <j>X + 2 (/ 3 , - 45ff 2 / 4 ) X 2 + 2H (l5/ 4 ,0 - UH 2 f 5 ) j>X 2 + A2H 2 f 5 ^X 3 = , (7) 



'pi 

E 2 = 3M 2 } H 2 + X-V + 2 (m p 2 j - 6f 4 X 2 - 4Hf 5 <j)X 2 + 2/ 5 , X 3 ) H 
-2 [/ 3 + - (5/ 4 ,0 - 10H 2 f 5 ) X - 6Hf 5 ^xj Xij> 

~2 (/ 3 ,0 + 9H 2 U) X 2 -A (3/ 4 ^ + 2H 2 h) H^X 2 + 2 (2/ 4 ,00 - H 2 f 5 ^) X 3 + AHf^X 3 = , 
E 3 = 3#0 + V^ + l8H 2 f 3 X + 108H 3 U<j)X - 2 {f iM + 18H 2 fa - 30if 4 / 5 ) X 2 

+ + 6Hf 3 <P - 4 (/ 3 ,0 - 27H 2 U) X - 20 (3/4,0 - 2H 2 f 5 ) H^X - 90H 2 / 5 . X 2] 



(8) 



3/3 + 36Hf 4 j> - (15/4,0 - 30F 2 / 5 ) X - 18Hf 5 .^X 



XH 



-2 (12/4 



19^ 2 /5,0j Hcf>X 2 - ZQH 2 f 5M X 



0. 



(9) 



where H = a/a is the Hubble parameter, and a dot denotes a derivative with respect to t. 

Let us consider the covariant Galileon Lagrangian where the functions fa (i = 3,4, 5) are given by Eq. In this 
case it is convenient to introduce the following dimensionlcss quantities 



pi 



MM 



pi 
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and 



Mi, 



r 



M 2 M 2 pl ' 



M*0 



M 2 M p i 



M 



The constraint equation ([7]) can be written as 

6z 2 -y 2 - 2U(x) - 6c 3 ay 3 z + 45c 4 a 2 y 4 z 2 - A2c 5 a 3 y 5 z 3 = . 
Combining Eqs. ([7]) and © to eliminate V and then using Eq. ([§]) to solve for <f> and H, it follows that 



da; 
dr 



dz 
oV 



(10) 



111) 



(12) 



(13) 



[9c 2 3 a 2 y b z + 3c 3 ay 2 {78c 5 a 3 y b z 3 - 63c 4 a 2 y 4 z 2 + y l - 6z 2 ) + 810c 2 4 a 4 y 7 z 3 - 3c 4 a 2 y 3 (603c 5 a 3 y 5 z 4 



+15y 2 z + U$(x)y - 36z 3 ) + 9A5c 2 5 a b y 9 z b + 3c 5 a 3 y 4 z(2ly 2 z - 30z 3 + 2Uj(x)y) - 2U, 4> {x) - 6yz}/A , (14) 



where 



= -[27c 2 3 a 2 y 2 z 2 + c 3 a{A50c 5 a 3 y 4 z 4 - 432c 4 a 2 y 3 z 3 + \2yz + Uj(x)) + I620c 2 4 a 4 y 4 z 4 - 12c 4 a 2 yz{9yz 
+270c 5 a 3 y 4 z 4 + U i<t> {x)) + 1 + 1575c 2 5 a b y e z e + lbc b a 3 y 2 z 2 {8yz + U^{x))]y 2 /A , 



2 + 3c 3 ay(4z + c 3 ay 3 ) - 3c 4 aV(36z 2 - y 2 + \8c 3 ay 3 z - 90c 4 a 2 y 4 z 2 ) 
+3c 5 a 3 y 3 z(40z 2 - 2y 2 + 18c 3 ay 3 z - 192c 4 a 2 y 4 z 2 + 105c 5 a 3 ?/ 5 z 3 ) . 



(15) 



(16) 



Numerically it is usually more stable to solve Eqs. (|T3"|) and (|T4"]) with the constraint equation (|T2"j) rather than solving 

Eqs. mg-cis). 



B. The spectra of density perturbations 



The spectra of scalar and tensor perturbations generated in the theories given by the action ([T]) were derived in 
Refs. |31rl33j . Here, we briefly review their formulas in order to apply them to concrete inflaton potentials. 
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The perturbed line element about the flat FLRW background is given by [4( 



ds 2 = -(1 + 2A)dt 2 + 2d i Bdtdx i + a 2 (t) [(1 + 2K)5 i: 



hij] dx l dx° 



(17) 



where A, B, 7Z are scalar metric perturbations, and are tensor perturbations which are transverse and traceless. 
The inflaton field is decomposed into the background and inhomogeneous parts, as (f) = 0o(O + S<p(t,x). We choose 
the uniform- field gauge characterized by 8<j) = 0, which fixes the time-component of a gauge-transformation vector 
The scalar perturbation E, which appears as the form Ejj in the last term of (fl7|) . is gauged away, so that the 
spatial part of £ M is fixed. Vector perturbations decay during inflation, so that their contribution is negligibly small. 

We expand the action ([T]) up to second-order in perturbations by using the Hamiltonian and momentum constraints. 
For the theories given by Eqs. (U|) and (|5|) the second-order action for scalar perturbations reduces to [U HH 



S {2) = / dtd 3 xa 3 Q s 



K 2 - ^{&Rf 



where 



and 



Q s = 



Wi(4wiW3 + 9w 2 ) 



3w 2 



i(2w\w2H — W 2 W4 + 4wiWiW2 — 2w\lD2) 
W\{4:WlW3 + Qw 2 ) 



U>2 

w 3 
U)4 



2M 2 ,H 



-9M^H 2 



Ki- 



rn 

(19) 
(20) 

UH 2 f 5 ^j X 2 + 28Hf 5 ^X 3 , (21) 
f 3 <f>) X + 6 (l35H 2 h - 2/3,^ - 45ff/ 4 > + 56i? 3 / 5 0) X 2 - 504H 2 f 5 ^X 3 , (22) 
2/5, <t>X 3 . 



wx = Mpj -2 3/4 + 2Hf 5 </> X z + 2f 5 ^X 



2h4>X - 2 30fT/ 4 - 5/4 



2 / 4 -2/ 5 X 



(23) 



The conditions for the avoidance of ghosts and Laplacian instabilities correspond to Q s > and c 2 > 0, respectively. 
The two-point correlation function of the curvatureperturbation 1Z can be derived by employing the standard method 
of quantizing the fields on a de Sitter background |4C| ■ The power spectrum of the curvature perturbation TZ is 



H 2 



^ 2 Q s c 3 



(24) 



which is evaluated at c s k — aH (where k is a comoving wavenumber). 

We decompose the intrinsic tensor perturbation hij into two independent polarization modes, as hi 
Then the second-order action for tensor perturbations is given by 



s t 2) = ^J dtd3xa3 Qt 



5 (£ 



where p 



x , and 



n - Wl 



W4 
Wl 



We require that Q t > and c 2 > to avoid ghosts and Laplacian instabilities. The tensor power spectrum is 



Vi 



H' 2 



2^Q t c 3 



(25) 



(26) 



(27) 



which is evaluated at c±k = aH. 



C. Slow-roll analysis 

For the covariant Galileon theory Q we employ the slow-roll approximation to estimate the physical quantities 
introduced in previous subsections. We define the following slow-roll parameters 

H X c^X 2c 4 X 2 _ 6c 5 Hj>X 2 
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whose magnitudes are much smaller than 1 during inflation. Eliminating the term V from Eqs. (JT)) and ([5]). we obtain 
the following equation 

e = : ^ ^ 6 X + 3d 3 + 18d 4 + 5d 5 . (29) 

1 — 604 — Z65 

In the second approximate equality we neglected the terms at second-order in slow-roll. 
Under the slow-roll approximation the field equations ([7]) and © reduce to 

3M^H 2 ~ V , (30) 
3fty(l + -A) + U ~ , (31) 



where 



, , Hj> ( H A\ A , (HA* 35 3 + 186 4 + 56 5 

A = 3c *M3- 1SCi [w +15C5 (w) = r x ■ (32) 



Using Eqs. (|3ll|) and (f3"Tj) . the parameter <5jf can be estimated as 



where 



From Eqs. ^SJ and (02) it follows that 



5x ~ (T+Ap ' < 33 > 



(1 + ^-^. (35) 



The conventional slow-roll inflation corresponds to the limit A — > 0, in which case e ~ ~ In the regime where 
is much larger than 1 the evolution of the field slows down relative to that in standard inflation. 
We define the number of e-foldings from the time t to the time tf at the end of inflation, as N = J. f H(t) dt. From 
Eqs. (J30J) and §5T§ we have 

N — t?2" / [l + A)—d$. (36) 
M P i 

The field value <pf at the end of inflation is known by solving e(</>/) = 1, that is 

e^f) = 1 + A(cf> f ). (37) 

Since the factor A in Eq. (|32|) involves the field velocity, we need to express 4> in terms of 4> according to Eq. (|3T|) for 
the evaluation of </>/ and N. 

Under the slow-roll approximation the quantities Q s and c 2 read 



Q s ~ Af£ (J* + 6^3 + 54<5 4 + 20<5 5 ) , (38) 

(39) 



n2 _ fa + 4(5 3 + 26(5 4 + 865 
Zs ~ 5 X + 663 + 5464 + 2065 



In the regime where the Galileon self- interactions dominate over the standard kinetic term we have {|^3|, \8&\, \6s\} ^> 
6x- In order to avoid that Q s becomes negative we demand the following conditions 

c 3 > , c 4 < , c 5 > . (40) 

If 6x is much larger than |<$3|, |<5 4 |, and |#5|, then the scalar propagation speed squared is close to 1. If either of 5, 
(z = 1,2,3) is the dominant contribution in Eq. (|39[) . we have 

c 2 ~ 2/3 (63 dominant) , (41) 

c 2 ~ 13/27 (<S 4 dominant) , (42) 

c 2 ~ 2/5 (J 5 dominant). (43) 
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This shows that the Laplacian instability of scalar perturbations is absent during slow-roll inflation. 

-t 



The quantities Q t and cf are approximately given by 



M 2 

Q t ~ -^(1- 3(5 4 -2<5 5 ), c 2 ~ f + 4^ + 2^5, (44) 

which are both positive. 

Under the slow-roll approximation the power spectra of scalar and tensor perturbations are given, respectively, by 

H 2 1 V (dx+6S 3 + 54S i + 205 5 ) 1 ^ 

3 ~ 8w 2 M^ c s e s ~ 24n a M* 1 + 4<J 3 + 26<5 4 + 8<5 5 ) 3/2 ' 1 j 
2_£/"^ 2V^ 

".-^^ < 46 > 



where 



The tensor-to-scalar ratio is 



00- ~ 5 X + 4(5 3 + 26<5 4 + 8<5 5 . (47) 
M P i 



Vt ir 1fi fe+ 4,5 3 + 26(54 + 8(5 5 ) 3 / 2 

r = n = 16Cs£s = 16 M6^TW^ • (48) 



Defining the spectral indices as n s — 1 = <fln'P s /cnnfc| Cs fc =a ij and n t = enn'P t /<ilnfc| Ct ; c=Q #, it follows that 



where e is given in Eq. ([29]) . and 



n s - 1 = -2e-r] s -s, (49) 
«t = -2e, (50) 



^ = ife? S= il/ (51) 
The consistency relation between r and nt is 

r = -8c s (n t - 25 3 - 16<5 4 - 6<5 5 ) ■ (52) 

In the regime Sx \6i\ (i = 1,2,3) the standard consistency relation r = —8n t holds. If cither of the terms \8i\ 
(i = 1, 2, 3) dominates over other terms, it follows that 

r = — 8.71nt ((5 3 dominant) , (53) 
r = -8.02n t (64 dominant) , (54) 
r = — 8.10n t (i5 5 dominant) . (55) 

Since the ratio r/n t is close to —8 in all cases, the observational bounds on n s and r are similar to those derived by 
using the consistency relation r = —8nt- 

III. THEORIES WITH G 3 / 0, G 4 = 0, G 5 = 

We first study the covariant Galileon theory in which only the term — (c 3 /M 3 )XD</> is present in the action (TTJ, i.e., 

c 3 ^0, c 4 = 0, c 5 =0. (56) 
Solving the slow-roll equation (|3ip for (j), it follows that 



1-^-1 , *nd AW = - . (57) 



6c 3 # IV M 3 J ' vv ^ 2 I V M 3 
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For c 3 > one has > and < from Eqs. (gOJ) and (|3T|). If c 3 < 0, then < and V<p > 0. In the former 
and latter cases we choose the coefficients c 3 = 1 and c 3 = —1, respectively, without loss of generality. The transition 
from Galileon inflation to standard inflation can be quantified by the condition A(4>g) = 1, which translates into 

c 3 V^ G ) = -2M 3 . (58) 

The field value </>/ at the end of inflation is known from Eq. (j3"7)) . i.e., 

1 



V AP 



(59) 



For the scalar potential with V$ > the transition from the regime 5 3 3> 8x to the regime (5 3 <C Sx occurs during 
inflation provided that \4>g\ > On the other hand, if \4>a\ < \<fif\, the Galileon self- interaction dominates over 

the standard kinetic term during the whole stage of inflation. 

Since c 3 V^/M 3 = —.4.(1 + .4), the number of e-foldings (|3"6")) reads 

Using Eqs. (|33|) . (f34|) . and the relation 5 3 = (.4/3)<5a', the scalar power spectrum (|45|) reduces to 

U 3 (l+.4) 2 (l + 2.4) 1 / 2 



12tt 2 M^V% (1 + 4.4/3) 3 / 2 



(61) 



For a given inflaton potential and a mass scale M , the field value 0/ is known by solving Eq. (|59p . Then the number 
of e-foldings can be evaluated from Eq. to find the value of (/) at ./V = 60 (for which we denote (/>6o)- The WMAP 
normalization of the scalar power spectrum is P s {4>6o) = 2.4 x 10~ 9 Q, by which the mass M can be related to model 
parameters of the potential. 

The scalar spectral index n s is known from Eq. (|6ip according to n s — 1 = T' s /{H'Ps)- Taking the time derivative 
of the quantity A given in Eq. ([57)1 and using Eq. (|3"0"j) . it follows that A/H = —rj^A/il + 2A), where 

V, = M^. (62) 



V 



Then we obtain 



6e0 27^ 



1 + A 1 + 4.4/3 



.4 



6(1 + 2^) : 



(63) 



where e<p is defined in Eq. (j34|) . For .4 — > the formula (|63|) recovers the result n s — 1 ~ — 6e^ + 2ri0 in conventional 
slow- roll inflation. In the limit A ^ 1 we have that n s — 1 ~ — 6e</,/.4 + 377^/(2^4). From Eq. (|48| the tensor-to-scalar 
ratio reads 

r = We, (^^Z 3 ) 372 (64 ) 

^ (1+^) 2 (1+ 2^)1/2- ^ 

In the limit A — > this reproduces the standard relation r ~ 16e0, but for A 3> 1 we have r ~ 64a/6 £0/(9.4). The 
observables (j6"3")l and are functions of 4>(N), so that they can be evaluated for a given inflaton potential. 

If the Galileon term is dominant even after the end of inflation, this affects the oscillation of inflaton during 
reheating. In order to avoid that the 1 + 6H / 3 </) term in front of (f> in Eq. §§§ become negative, we require 

1 + 6Hc 3 0/M 3 > . (65) 



The field velocity <j) changes its sign during the oscillating stage of inflaton. This means that the condition (p5|) 
can be violated depending on the model parameters. Note that the determinant (|16[) is approximately given by 
A ~ 2(1 + 6Hc3<p/M 3 ), so that the violation of the condition (|65|) leads to the divergence of Eqs. (fT4|) and (fTB"]) 1 . In 



1 Note that a similar determinant singularity appears in the context of anisotropic string cosmology [4 if . 
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the regime 1 + 6Hc3<j)/M 3 < the field climbs up the potential like a phantom field, so successful reheating cannot 
be realized. 

For smaller values of M the condition (|65[) tends to be violated. In this case we also show that c 2 s can be negative 
due to the dominance of the Galileon term. For the potentials (a) V{4>) = X(f> n /n (chaotic inflation) and (b) V(4>) = 
A 4 [l + cos((/)//)] (natural inflation), we clarify the parameter space in which the oscillation of inflaton is not violated 
and c 2 s remains positive. We also place observational bounds on each model from the information of the scalar spectral 
index n s and the tcnsor-to-scalar ratio r. 



A. Chaotic inflation 



First we study the case of chaotic inflation characterized by the potential 

V&) = V , (66) 
n 

where n and A are positive constants. The initial value of the field (j> is assumed to be positive, so that <fi < and 
hence C3 = —1. From Eq. ([58]) the field value at the transition is given by 

G =(2M 3 /A) 1/( "- 1) . (67) 

If the slow- roll parameter (|35[) at <f> = <f>Q is smaller than 1, the transition from Galileon inflation to standard inflation 
occurs during inflation. The condition c(</>g) < 1 translates into 

M > 2-™/V^ 1 )/ 3 M^- 1)/3 A 1 / 3 = M c . (68) 

For the potential ([66]) the function A(<fi) in Eq. (|57|) is given by 



in which case, apart from the case n = 2, the number of e-foldings (|60|) is not integrated analytically. Moreover, 
in order to find the field value </>/, we need to solve Eq. (|59|) numerically. In the limit A 3> I, however, it is 
possible to derive the analytic expression of <j> in terms of N. Since A((f>) ~ -J A0™ _1 /A/ 3 in this limit, we have 
l Ml A/ 3 / 2 /(2v / A) from Eq. ([59]). Then, integration of Eq. ([60]) gives 



6 (n+3)/a „ 



, J _„ W0 nM^M 3 / 2 
^(n+3)/2 „ pi — [(n + 3)A r + n j _ (70) 

2y A 

Substituting this solution into Eq. (|61[) and using the WMAP normalization T s = 2.4 x lfr 9 at N = 60, we find 

M p 4 ! " (61n + 180)"+! ' 1 j 

For smaller M, A tends to be larger. The scalar spectral index ([63]) and the tensor-to-scalar ratio ([64]) reduce to 

(n + 3)Ar + n' 9 (n + 3)7V + n ' y ' 

which agree with those given in Ref. [23|. For n — 4 and N = 60, for example, we have n s = 0.965 and r = 0.164. 
In the limit i«lwc have cj> } = nM pl /y/2 and 2 = 2nM 2 l (N + n/4). The WMAP normalization at = 60 gives 

A = 2.8 x 10- 7 n 3 [n(120 + n/2)]^ ( " +2)/2 M 4 f " . (73) 

The scalar spectral index and the tensor-to-scalar ratio are 

2(n + 2) 16n 
n s = 1 v y , r = , (74) 

which correspond to those for standard chaotic inflation. In the regime between A ^ 1 and A <C I we need to 
evaluate n s and r numerically For n — 4 and A = 60, for example, we have n s = 0.951 and r = 0.262. 
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Figure 1: The mass parameters M and m satisfying the WMAP normalization V a = 2.4 x 10 at N = 60 for the quadratic 



potential V((f>) = rn 2 (f> 2 /2 with the term G3 
not violated during reheating. 



-X/M . The solid line represents the region in which the inflaton oscillation is 



V(ct>) = m 2 2 /2 



Let us consider the case of the quadratic potential V{4>) = m 2 (j) 2 /2, i.e., A 
WMAP normalization (f7T|) leads to the following relation 



m 



M 



M p i M p i 



4.1 x 10" 



and n = 2. In the limit A 3> 1 the 



(75) 



In another limit ^4 1, we have m ~ 6.2 x 10 -6 M p i from Eq. ([75|. In Fig. Q] we plot m versus M constrained by the 
normalization V s = 2.4 x 10~ 9 at N = 60. For small M satisfying M/M p \ <C 10 -3 the numerical result in Fig. [1] is in 
good agreement with the analytic estimation (|75|) . In the regime M/M p \ 3> 10~ 3 the mass m approaches the constant 

value m ~ 6.2 x 10~ 6 M p i. Under the condition ([55]). i.e., M > 2 _1 / 3 M p L 1 /3 TO 2 / 3 , the transition from the regime .4 > 1 
to the regime A < 1 occurs during inflation. Combining this condition with the constraints on M and m shown in 
Fig- H it follows that M > 4.0 x 10 _4 M p i. If M < 4.0 x 10 _4 il/ p i, the Galilcon self-interaction dominates over the 
standard kinetic term during inflation. 

In order to see the effect of the Galilcon term during inflation and reheating, we numerically solve the background 
equations (fT2"]) - (|14[) with the initial conditions determined by the slow-roll analysis. We confirm that the slow-roll 
approximation is accurate enough to reproduce the numerical values of N with the difference less than a few percent. 

In Fig. [5] we plot the evolution of <j> and c 2 s for three different mass parameters M and m constrained by the WMAP 
normalization. The case (a) corresponds to the mass M = 3.0 x 10 _4 M p i, which is smaller than the critical mass 
M c = 4.0 x 10 _4 Af p i. Hence the Galileon self-interaction dominates over the standard kinetic term by the end of 
inflation. As we see in the right panel of Fig. [2 the solutions enter the regime in which c 2 s is negative. For M smaller 
than 3.0 x 10 _4 M p i the period in which c 2 is negative tends to be longer with \c 2 \ much larger than 1. Since scalar 
perturbations grow very rapidly in such cases, the Universe becomes inhomogeneous at the level of destroying the 
homogenous background. The case (b) shown in Fig. [5] corresponds to the marginal one in which the minimum value 
of c 2 is 0. In the case (c) the transition from the regime A > 1 to the regime A < 1 occurs during inflation and c 2 
always remains positive. For the range of masses M used in the numerical simulations of Fig.[5J the inflaton oscillation 
is not violated as long as the backreaction of created particles is neglected. 

The condition for the avoidance of negative values of c 2 is 



M > 4.2 x 10 _4 M 



P l 



(76) 
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Figure 2: Evolution of the field <f> (left) and the scalar propagation speed squared c 2 (right) for the quadratic potential 
V(<j>) = m 2 4> 2 /2 with the term G 3 = -X/M 3 in three different cases: (a) M = 3.0 x 10~ 4 M p i, m = 1.45 x 10~ s M p i, (b) 
M = 4.2 x 10~ 4 M p i, m = 1.1 x 10" 6 M p i, and (c) M = 1.0 x 10~ 3 M p i, m = 6.9 x 10" 6 M p i. We choose the initial conditions 
at A*' = 60 determined by the slow-roll analysis, i.e., (a) Xi = 6.28, yt = —5.25 x 10~' ! , z% = 1.24 x 10 _1 , (b) Xi = 8.32, 
yi = -5.07 x 10 -3 , Zi = 8.89 x 10~ 2 , (c) x t = 13.55, y { = -3.90 x 10 -3 , Zi = 3.82 x 10~ 2 , respectively. In the case (a) the 
system enters the region with negative values of c? s , whereas in the case (c) (? s is always positive. The case (b) is the marginal 
one in which the minimum value of (? a is 0. 
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Figure 3: Observational constraints on the quadratic potential V((f>) = m 2 (f> 2 /2 with the term G3 = —X/M 3 in the (n s ,r) 
plane for the numbers of e-foldings N = 50 , 60 , 70. The thin solid curves show the lcr (inside) and 2<r (outside) observational 
contours constrained by the joint data analysis of WMAP7, BAO, and HST. For smaller values of M the tensor-to-scalar ratio 
r gets smaller, whereas the scalar spectral index n a increases. 
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under which c 2 . finally approaches 1 without entering the regime c 2 < 0. Note that the conditions c 2 > 0, Q s > 0, 
and Q t > are satisfied for the above values of M. Numerically we find that the inflaton oscillation is not violated 
during reheating for 

M > 2.5 x l(T 4 A/ pl , (77) 

which is related to the condition (|65[) . During inflation in which c^cj) is always positive, the condition (|65[) is always 
satisfied. However, after <f> changes its sign during reheating, the condition ([65)) is violated for M < 2.5 x 10 -4 M p i. 
The criterions ([76")) and ([77)) are not very different from each other. 

In Fig. [3] the theoretical values of n s and r are plotted as a function of M ranging in the region (|77|) with three 
different values of N (= 50, 60, 70). We also show the la and 2a observational contours constrained by the joint data 
analysis of WMAP7 pj, Baryon Acoustic Oscillations (BAO) [42[, and the Hubble constant measurement using the 
the Hubble Space Telescope (HST) [43j . As we decrease the value of M, the two observables shift from the values 
in Eq. f|74|> to those in Eq. (|72p. For smaller M, r gets smaller whereas n s increases, so that the quadratic potential 
shows better compatibility with the data. Even for the mass M corresponding to the lower limit of Eq. ([77)1 the term 
A is larger than 1 during most stage of inflation, in which case n s and r are close to the asymptotic values given in 
Eq. (dl). 



2. V(<t>) = \<j) 4 /4 

We proceed to the case of the quartic potential V{4>) = \<f) A /4. In the limits A ^s> 1 and i « 1 the WMAP 
normalizations (|7T|) and ([73)) give \(M/M p i) 4 ~ 2.4 x 10 -26 and A ~ 1.6 x 10~ 13 , respectively. Figure 4 shows the 
viable parameter space in the (M, A) plane satisfying the WMAP normalization at N = 60. Under the condition 
([55]). i.e., M > 2 2 / 3 A 1 / 3 Af p i, the transition from the regime A > 1 to the regime A < 1 occurs during inflation. 
Combining this condition with the constraints on M and A shown in Fig. SJ it follows that M > 2.8 x 10 _4 M p i. If 
M < 2.8 x 10~ 4 M p i, the Galileon term dominates over the standard kinetic term during inflation. 

By solving the background equations of motion (|12[) - ([14[) . we find that c 2 s remains positive for 

M > 1.7 x 10~ 4 M p i . (78) 

The inflaton oscillation is not violated during reheating provided that 

M > 9.5 x 10~ 5 M pl . (79) 

In Fig. [5] the theoretical values of n s and r are plotted as a function of M ranging in the region ([79)) with N — 
50, 60, 70. In the limit M — > oo the quartic potential is outside the 2a observational contour for N smaller than 70. 
In the presence of the Galileon term the model can be compatible with the current observations due to the suppressed 
tensor-to-scalar ratio and the larger scalar spectral index. For N = 60 the model is within the 2a contour for 

M < 7.7 x 10~ 4 M p i . (80) 

In terms of the parameter A the conditions ([79)) and ([80)) translate into 

3.4 x 10~ 13 < A < 3.1 x 10~ 10 . (81) 

Under the constraint ([78)) the upper bound is A < 3.0 x 10~ n . The result ([8T)) shows that one cannot accommodate 



the self coupling A ~ 0.1 of the Higgs boson in the presence of the coupling G3 = —X/M 3 . 
We also studied the case of the generalized Galileon term —6*3(0, X)D<j), where 

G3 = ^X , (82) 

which was proposed in Ref. [27| ■ Numerically we find that the violation of the inflaton oscillation does not occur for 
M > 3.6 x 10" 4 Af p i and A < 2.7 x 10" 8 . The self coupling A is still much smaller than the order of 0.1. 

B. Natural inflation 

Natural inflation [2(| is characterized by the potential 



V{<j>) = A 4 



(83) 
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Figure 4: The parameters M and A satisfying the WMAP normalization V s = 2.4 x 10 at TV = 60 for the quartic potential 
V((f>) = \(f> 4 /4 with the term G3 = —X/M 3 . The solid line represents the region in which the inflaton oscillation is not violated 
during reheating and the model is within the 2a observational contour in the (n a ,r) plane. 




Figure 5: Observational constraints on the quartic potential V{<j>) = A</) 4 /4 with the term G3 = —X/M 3 for three different 
values of TV. The la and 2a observational contours are the same as those in Fig. [3] While the standard case (M — > 00) is 
outside the 2a bound, the presence of the Galileon term can make the quartic potential compatible with observations. 



where A and / are constants having the dimension of mass. In the absence of the Galileon term this potential can be 
compatible with the observational data only for / > 3.5M p i [2lj . This is the regime in which standard quantum field 
theory may not be reliable. If the field <fi is a string axion, / is usually smaller than the order of M p i [22|, 13 ■ In the 
following we study whether this problem can be alleviated or not in the presence of the Galileon term G3 = C3X/M 3 . 
We assume that inflation occurs in the region < <p/ f < n, in which case </> > 0. We choose C3 = 1 to satisfy the 
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Figure 6: Observational constraints on natural inflation with / = 0.1M p i in the presence of the term G3 = X/M 3 for three 
different values of N. The parameter range of 7 = A 4 /(M 3 A'/ P i) corresponds to 4.0 x 10 5 < 7 < 1.0 x 10 11 . The la and 2a 
observational contours are the same as those in Fig. [3] 



condition (|40p. For the potential (|83p the function A((f>) in Eq. ([57)) is given by 



1 / / 4 7 sin(0//) 



■4W = -h/l+ ' J - 1 h (84) 



where 



9= i4' 7= ^v (85) 

Note that in the limit (j> — > one has A — > 0. We are mainly interested in the case where the initial displacement of 
the field 4>i satisfies the condition A{4>i) > 1. This can be achieved for 47 3> <j provided that ^ is not very close to 0. 
The field value </>g at the transition (.4. = 1) is given by 

sm(<f> G /f) = 2q/~f. (86) 

For the existence of 4>g we require that 2q < 7, i.e., M 3 < A 4 /(2/). From Eq. (JSHJ) the field value </>/ at the end of 
inflation satisfies 



1 - cos^///) = 2 
1 + cosfa///) 




i 47sin(0£//) 
9 



(87) 



The transition from Galileon inflation to standard inflation occurs under the condition <f>Q <</>/. In the limit 7—5-0 
one has cos(</>///) = (1 — 2q 2 ) / (1 + 2q 2 ) , so that <f>f/f — >• for q -C 1. This implies that, in the absence of the Galileon 
self-interaction, it is difficult to realize sufficient amount of inflation for / <C M p \. If the Galileon term is present with 
7 — > 00, the field value </>/ can be close to nf even for / <C M v \. 
The slow-roll parameters and 77^ are given by 

1 sin 2 (0//) 1 cos(0//) 

* 2g 2 [l + cos(0//)] 2 ' V * q 2 l + cos(0//) ' 1 j 

In the limit 7—5-0, i.e., A — >• 0, the scalar spectral index is n s ~ 1 — 6e^ + 277^. When q <C 1 inflation occurs in the 
region 0// <C 1, so that \r)$\ ~ 1 /(2g 2 ) 3> 1. This case is in contradiction with observations because n s significantly 
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Figure 7: The allowed values of 7 versus f/M p i in natural inflation with the term Gs = X/M 3 . In the region (i) the model is 
within the 2a observational contour in the (n s ,r) plane for iV = 60. In the region (ii) the inflaton oscillation is not violated 
during reheating. There is a viable parameter space only for f/M p i > 1.7. 



deviates from 1. In another limit 7 — > 00 one has A — > 00 and hence the scalar spectral index (|63|) can be as close as 
1 even for g« 1. In this limit, inflation occurs in the regime close to the potential minimum ((f)/ f = tt). Since the 
potential is approximately given by V(<p) ~ (A 4 /2/ 2 )(0 — irf) 2 in this regime, n s and r arc the same as those given 
in Eq. ([72]) with n = 2. Hence, for 7 — > 00, it follows that 

9 128V6 
n s = 1 , r = — — - , 89 

bN + 2 1 9(5iV + 2)' V ' 

which give n s = 0.970 and r = 0.115 for N = 60. 

In the intermediate regime characterized by < 7 < 00 we need to evaluate n s and r numerically according to 
Eqs. (f6"3")) and For given values of /, A, and M, cj)f is known by solving Eq. (|87[). Integrating Eq. (joT))) numerically, 
we can determine the field <fi in terms of the number of e-foldings N. The WMAP normalization V s = 2.4 x 10~ 9 at 
N = 60 provides one constraint between the three parameters /, A, and M . In other words, for a given /, the two 
parameters A and M are related to each other. 

Let us first consider the case / = 0.1M p i, i.e., q = 0.1. In Fig. [51 we plot the theoretical values of n s and r in 
the range 4.0 x 10 5 < 7 < 1.0 x 10 11 for three different values of N, together with the la and 2a observational 
contours. When N = 60 the model is within the 2a contour provided that 7 > 4.3 x 10 5 . In the limit 7 — > 00, n s 
and r approach the asymptotic values given in Eq. (|8"9")l . This asymptotic case is within the la contour for N > 55. 
For larger 7, however, the inflaton oscillation during reheating tends to be disturbed by the Galileon term. After the 
inflaton velocity </> changes its sign from positive to negative, the condition (|65|) can be violated for 7> 1. Detailed 
numerical simulations show that the violation of inflaton oscillations does not occur provided that 7 < 0.5. Thus 
there is no viable parameter space of 7 satisfying both the WMAP bound and the successful reheating. We also note 
that, if 7 is larger than 0.05, c 2 s becomes negative. Hence, for 7 > 0.5, the model is plagued by the reheating problem 
as well as the negative instability of scalar perturbations. 

We also study the models of other values of / ranging in the region 0.1 < f/M p \ < 2. In Fig. [7J we show the 
two kinds of constraints on the parameter 7 versus f/M p \. Above the dotted line (i) the model is within the 2a 
observational contour in the (n s ,r) plane, whereas under the solid line (ii) the oscillation of inflaton is not violated 
during reheating. For the compatibility of two constraints we require that / is bounded to be 

/>1.7Mpi. (90) 



Hence the problem of the super-Planckian values of / in standard natural inflation is not circumvented by the Galileon 
term G 3 = X/M 3 . 
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IV. THEORIES WITH G 3 = 0, G 4 + 0, G 5 = 



We proceed to the covariant Galileon theory (|6]) with 

C3 = , c 4 / . 



c 5 = 0. 



(91) 



Since C4 < to avoid ghosts, we set C4 = — 1 without loss of generality. From Eq. (j3"Tj) the quantity A = l8(H<j>/Al 3 ) 2 
satisfies the following relation 

A(l+A) 2 = 2V 2 4> /M 6 . (92) 
The field value <pQ at the transition from Galileon inflation to standard inflation obeys 

V%{<fa) = 2M 6 . (93) 
From Eq. (|3"2"]l we have 64 = ASx/^8. Using Eq. (|33[) . the scalar power spectrum (|4"5]) can be written as 

V 3 (1 + 4) 2 (1 + 34)^2 



8 127r 2 Ai£y| (i + i3^/9) 3 / 2 ' 

Taking the time derivative of Eq. ([92]) and making use of Eq. (|3"Tj) . it follows that A/H 
scalar spectral index is given by 



71, - 1 = - 



The tcnsor-to-scalar ratio (l48l) reads 



6e^ 2??0 
1+4 1 + 3.4/2 



1 - 



34(5 + 8A) 



2(l + 34) 2 (9+134) 



166. 



(1 + 13.4/9) 



3/2 



(1+4) 2 (1+ 34)^2 



(94) 

277^/(1 + 3-4). Then the 
(95) 

(96) 



In the limit i->oowe have n s - 1 ~ -66^/4 + 4770/(34) and r ~ 208 V39 6^/(814). 

In the following we focus on the potential V(cf>) — Xcj) n /n of chaotic inflation. From Eq. the field value at the 
transition is 

G =(2M 6 /A 2 ) 1/[2( - 1)] . (97) 
The condition under which the transition occurs during inflation corresponds to £(4>g) < lj which translates into 

(98) 



M > 2 ( 1 -2»)/6 n («~l)/3 M ^-l)/3 A l/3 



Let us consider the case in which the condition A ^S> 1 is satisfied during the whole stage of inflation. Then we 
have 54(77 <j>) 3 ~ —M 6 V^ from Eq. (|31[) . The end of inflation is characterized by the condition e(cj)f) ~ e ( p((f>f)/A = 1, 
which gives 4>f — [n M 3 l M 3 /(4X)} 1 ^ n+2 \ The number of e-foldings (J36j) is related to the field value $ during inflation, 
as 



^ (n+2)/3 -^S[4(" + 2)^ + 3n]. 
The WMAP normalization V s = 2.4 x 10~ 9 at N = 60 provides the relation 

M 3n 2 3 "+ 2 71 n /2+ 4 (1.8 X 10~ 6 )"+2 



\ 2 



M 



n+8 
Pi 



39™/2+! (81n+ 160)( 5 ™+ 4 )/ 2 



Substituting Eq. ([99]) into Eqs. ((95j) and ([96]) in the regime 4 > 1, it follows that 

2(5n + 4) 



208V39 



4{n + 2)7V + 3n ' 



27 4(n + 2)iV + 3n 



(99) 



(100) 



(101) 



For n = 4 and TV = 60, for example, n s = 0.967 and r = 0.133. The tensor-to-scalar ratio is smaller than that studied 
in Sec. IIIII in the regime 4 1. 
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Figure 8: The parameter space satisfying the normalization V B = 2.4 x 10 -9 at N = 60 for the quadratic potential V(4>) = 
m 2 <f) 2 /2 (left) and for the quartic potential V(</>) = A0 4 /4 (right) in the presence of the term Gi = X 2 /M 6 . The solid lines 
correspond to the regions in which c 2 is positive and the observational constraints of n s and r are satisfied. 




In another limit A <C 1, we have the same relations as those given in Eqs. ([75)) and ([71)1 . In the intermediate regime 
between .4 > 1 and A -C 1 we need to solve the background equations (|12 p -(|14 p numerically in order to find the 
values of n s and r as well as the relation between A and M from the WMAP normalization. In Fig. [8] we show the 
parameter space for the two potentials V(<p) = m 2 </> 2 /2 and V(<p) = A</> 4 /4 satisfying the WMAP normalization at 
N = 60. In the two asymptotic regimes A 3> 1 and A <C 1, the analytic estimation given above agrees well with the 
numerical results. 

Unlike the case of the coupling G3 = C3X/M 3 , the term 1 + 54H 2 <j) 2 /M 6 in front of <fi in Eq. ((9]) remains positive 
even if <j> changes its sign. Numerically we confirmed that the determinant A defined in Eq. (fll))) does not cross even 



17 




Figure 10: Observational constraints in the (n s ,r) plane with the numbers of e-foldings N = 50,60,70 for the quadratic 
potential V((f>) = m 2 (f> 2 /2 (left) and for the quartic potential V(4>) = A0 4 /4 (right) in the presence of the term G 4 = X 2 /M 6 . 
The la and 2a observational contours are the same as those shown in Fig. [3] 



for the mass M much smaller than the r.h.s. of Eq. (|98|) . In the left panel of Fig. [9] we show the field evolution during 
reheating for the quartic potential V(<fr) = A</> 4 /4 with A = 0.1 and M = 7.3 x 10~ 6 M p i. In fact the the inflaton 
oscillation is not violated by the dominance of the term G4 = X 2 /M 6 . In this case, however, the scalar propagation 
speed squared oscillates significantly between largely negative and positive values (see the right panel of Fig. [5]). This 
leads to the strong enhancement of scalar perturbations for the modes inside the Hubble radius during reheating. 
While this instability does not directly affect the evolution of large-scale density perturbations relevant to CMB, the 
rapid growth of perturbations can invalidate the analysis without the backreaction of created particles after some 
stage of reheating [45j . As we see in Fig. [3J our numerical simulations without the backreaction effect show that c 2 
finally approaches 1 with oscillations. It remains to see how the created particles can change the evolution of <j> and 
c 2 at the late stage of reheating. This is beyond the scope of our paper, since nonlinear lattice simulations (along the 
line of Refs. |46[) are required to deal with such a problem properly. 

For larger values of M, the instability associated with negative c 2 tends to be less significant. For the power-law 
potential V(4>) = \(f> n /n we find that c 2 remains positive for 

M > 4.3 x l(T 4 M pl (for n = 2) , (102) 
M > 2.3 x 10~ 4 M p i (for n = 4) , (103) 

respectively. The regions in which these conditions are satisfied are shown as solid curves in Fig. |8] In Fig. [TO] we plot 
the theoretical values of n s and r for n = 2 and n = 4 as a function of M. Even for the lower bounds of Eqs. (|102p 
and (|103[) . n s and r are close to the values (|101[) corresponding to the limit „4 ^ 1. For the quadratic potential the 
presence of the term G4 = X 2 /M 6 leads to better compatibility with the WMAP data (see the left panel of Fig. IT0|) . 
In the case of the quartic potential the model is within the 2a observational contour under the condition 

M < 1.1 x 10~ 3 Af p i , (104) 

for N = 60. From Fig. [5] this condition translates into 

A>1.7xl0~ 13 . (105) 

If we demand the condition (|103j) for the avoidance of negative values of c 2 , the self coupling is bounded to be 
A < 9.9 x 10 _u . Recall that we do not have a constraint coming from the violation of inflaton oscillations. 
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V. THEORIES WITH G 3 = 0, G 4 = 0, G 5 + 

Finally we study the covariant Galileon theory ((6]) with 

c 3 = 0, c 4 = 0, c 5 ^0. (106) 
From Eq. ((HJ) the quantity A = 15c 5 i? 3 3 /M 9 satisfies 

A(l + A) 3 = -5c 5 Vl/(9M 9 ) . (107) 
The field value 4>g at the transition from Galileon inflation to standard inflation is determined by 

V%{<f> G ) = -72Af 9 /(5c 5 ) . (108) 
Since (S5 = ASx/5, the scalar power spectrum (|45|) reduces to 

p - y3 (i + -4) 2 (i + W /2 

/s 12^1^2 (i + 8^/5)3/ 2 ' { > 

On using the relation A/H = —3^^4/(1 + 4.4), the scalar spectral index is expressed as 



n. - 1 - 



1 + .4 1 + 8.4/5 
The tensor-to-scalar ratio (14"51) reads 



1 9A 



5(1 + 4.4) 5 



(110) 



(1 + 8.4/S) 3 / 2 

r=16 ^ (l + .A)2(l + 4^/ 2 - ( U1 ) 

Let us focus on the power-law potential V((j>) — Acff 1 /n. We assume that inflation occurs in the regime cf> > with 
the coefficient C5 — — 1 (under which the condition c§4> > is satisfied). Then the field value at the transition is 

G = [72M 9 /(5A 3 )] 1/[3(n - 1)] . (112) 

The condition under which the transition occurs during inflation is 

M > 0.94 • 2-"/V"- 1 )/ 3 A^r 1)/3 A 1 / 3 . (113) 

If «4 3> 1 during the whole stage of inflation, the field value at the end of inflation can be estimated as 4>f = 
[(9/5) 1 / 4 n 2 M2 1 M 9 / 4 /(2A 3 / 4 )] 4 /( 3 ™+ 5 ). The field <j> is related to the number of e-foldings N, as 

0( 3 n+ 5)/4 ^^$L( 9 -M1) [(3n + 5)N + 2n] . (114) 



4A 3 / 4 

From the WMAP normalization V s = 2.4 x 10~ 9 at N = 60 it follows that 

M$n „2(ra+5) 

A 5 = 1.33 x 10~ 7 ■ (4.44 x 10~ 7 ) Tl/2 • (1.36 x 10- 5 ) 3 "+ 5 — -7^-^ . (115) 

M 4("+ 5 ) v ; v ; (91n + 150) 7 "+ 5 v 

From Eqs. (|1 10[) and (|llip the asymptotic values of n s and r in the regime ^4^1 are 

7n + 5 256^10 n 
n s = 1 — , r = ; ; . 116 

{3n + 5)N + 2n' 25 (3n + 5)N + 2n v ' 

If n — 4 and JV = 60, for example, n s = 0.968 and r = 0.126. The tcnsor-to-scalar ratio is slightly smaller than 
that for the coupling G4 = X 2 /M 6 . In the regime « 1 the relations ([73")) and ([71)1 also hold for the coupling 
G 5 = -3X 2 /Af 9 . In the intermediate regime between A 3> 1 and ^4< 1 we resort to the numerical analysis to derive 
the relation between M and A from the WMAP normalization as well as to evaluate the observables n, and r. 
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In order to avoid that the term in front of <p in Eq. (O becomes negative, we require that 

1 + 120H 3 Xc 5 <j)/AI 9 > . (117) 

For smaller M this condition can be violated during reheating because of the sign change of <p (as it happens for 
the coupling G 3 = c 3 X/M 3 ). Note that this also leads to the divergence of Eqs. ([T4")l and (fT5")l through the crossing 
at A = 0. For the quadratic and quartic potentials we find that the infiaton oscillations are not violated under the 
conditions 

M > 2.7 x 1CT 4 7\/ P i (for n = 2) , (118) 
M > 1.5 x l(T 4 Mpi (for n = 4) , (119) 

respectively. The instability associated with negative values of c 2 tends to be stronger for smaller M, The conditions 
under which c 2 remains positive are given by 

M > 4.0 x l(n 4 A/pi (for n=2), (120) 
M > 2.9 x 10~ 4 M p i (for n = 4) , (121) 

respectively. The similar lower bounds on M (> 10 _4 M p i) also follow from Eq. (|113p . 

In the presence of the coupling G5 = — 3X 2 /M 9 the tensor-to-scalar ratio gets smaller relative to that in standard 
inflation, so that the quadratic potential V{<p) = m 2 (f> 2 /2 is compatible with the current observational data. For 
the quartic potential V{<f>) = X(p 4 /4: the model is within the 2a observational contour in the (n s ,r) plane under the 
condition 

M < 8.6 x 10~ 4 M pl , (122) 
for N — 60. Translating the conditions (|119l) and (I122[) in terms of the parameter A, it follows that 

1.6 x 10" 13 < A < 2.6 x 10 -9 . (123) 
As in the case of the coupling G3 = —X/M 3 , the self coupling A is required to be very much smaller than unity. 



VI. CONCLUSIONS 



We have studied the viability of potential-driven Galileon inflation described by the action (|T|). We mainly focused 
on the covariant Galileon theory in which the functions Gi (i = 3, 4, 5) are given by Eq. ([S]) with the choice ©. The 
Galileon self-interactions generally lead to the slow down for the evolution of the field, which allows the possibility 
to accommodate steep infiaton potentials. In Ref. [27j, for example, it was suggested that even the Higgs potential 
V(4>) = A</> 4 /4 with A ~ 0.1 can be consistent with the observed CMB temperature anisotropies because of the 
presence of the term G3 = C3X/M 3 . 

The dominance of the Galileon self-interactions relative to the standard kinetic term X can modify the dynamics 
of reheating after inflation. In order to clarify this issue, we numerically solved the background equations (|12l) - (fT5")) 
for several different infiaton potentials. We found that the oscillation of infiaton can be violated depending on the 
couplings Gi (i = 3,4,5). Moreover the dominance of the Galileon terms generally gives rise to the negative scalar 
propagation speed squared c 2 during reheating, which leads to the instability of small-scale density perturbations. 

For the theories where the covariant Galileon term G3 = C3X/M 3 is present, we found that the infiaton oscillation 
can be violated after the field velocity <j) changes its sign around the onset of reheating. This corresponds to the 
violation of the condition (p5"|) . which is related to the crossing of the determinant A in Eq. (|T6"|) at 0. The latter 
leads to the divergence of the background equations (|T4"|) and (|15[) . For the potentials V(<p) = \(f> n /n the violation of 
infiaton oscillations does not occur for M > 2.5 x 10~ 4 M p i (n = 2) and M > 9.5 x 10~ 5 M p i (n = 4). When n = 4 
this constraint translates into A < 3.1 x 10 -10 , which is much smaller than the coupling constant A ~ 0.1 of the Higgs 
boson. In the presence of the term G3 = C3X/M 3 the quartic potential V((f>) = A0 4 /4 is within the 2a observational 
contour in the (n s ,r) plane for M < 7.7 x 10~ 4 M p i. Taking into account this constraint, the self coupling is bounded 
to be 3.4 x 10~ 13 < A < 3.1 x 10~ 10 . We also found that c 2 s remains positive under the condition M > 1.7 x 10~ 4 M pl , 
which provides even the stronger upper bound A < 3.0 x 10~ n . We extended our analysis to the generalized Galileon 
term G3 = c^X/M 4 with the potential V(4>) = A</> 4 /4 and derived the bound A < 2.7 x 10~ 8 for successful reheating. 

In the presence of the term G3 = C3X/M 3 we studied the case of natural inflation described by the potential 
V(4>) = A 4 [l + cos(0//)] as well. While this potential can be compatible with the observed CMB anisotropies for 
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7 = A 4 /(M 3 M p i) ^> 1 even in the regime / <C M p \, the inflaton oscillation is violated under the condition 7 ^> 1. 
For the compatibility of two constraints, we found that / needs to be larger than 1.7 M p \. Hence the super-Planckian 
problem of the symmetry breaking scale in standard inflation (/ > 3.5M p i) is not improved significantly. 

For the Galileon coupling G4 = — c^X 2 /M & the scalar ghost is absent for C4 < 0, in which case the sign change of 
the determinant A in Eq. (|16[) can be avoided. In fact, we numerically confirmed that the oscillation of inflaton is not 
violated even for small M corresponding to the large self coupling A ~ 0.1 of the quartic potential V(<f>) = A</> 4 /4. On 
the other hand, for such small values of M, the scalar propagation speed squared c 2 heavily oscillates between largely 
negative and positive values (see the right panel of Fig. [5]) . This leads to the rapid growth of scalar perturbations for 
the modes inside the Hubble radius during reheating, which can invalidate the analysis without taking into account 
the backreaction of created particles. For the potentials V{<f>) = \<fi n /n the conditions for the avoidance of this 
negative instability are given by M > 4.3 x 10~ 4 M pl (n = 2) and M > 2.3 x 10" 4 A/ p i (n = 4). Taking into 
account this condition, the quartic potential V((j>) = A0 4 /4 is compatible with the current CMB observations for 
1.7 x 10~ 13 < A < 9.9 x KT U . 

In the case of the Galileon coupling G5 = 3c^X 2 /M 9 the sign change of A can occur for small M, as it happens 
for the coupling G3 = C3X/M 3 . For the potentials V(<fi) — \<fi n /n the inflaton oscillations are not violated for 
M > 2.7 x 10~ 4 M pl (n = 2) and M > 1.5 x 10~ 4 M pl (n = 4). Using the latter bound, the quartic potential is 
consistent with the CMB observations for 1.6 x 10~ 13 < A < 2.6 x 10~ 9 . We also found that the instability associated 
with negative c 2 is present for small M, which puts even severer upper bounds on A. 

Compared to the models of non- minimal field derivative couplings to the Einstein tensor [ljj [Tt}, the allowed 
parameter space of potential-driven Galileon inflation is more severely constrained because of the modified dynamics of 
reheating. We note, however, that there are some viable parameter spaces even for the quartic potential V((f>) = A(/> 4 /4 
due to the presence of the Galileon terms. It will be of interest to see whether future observations such as PLANCK 
[47j can place tighter constraints on such inflationary scenarios. 
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